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ABSTRACT 

We present a detailed study of the effects of gravitational microlensing on compact 
and distant 7-ray blazars. These objects have 7-ray emitting regions which are small 
enough as to be affected by microlensing effects produced by stars lying in intermedi- 
ate galaxies. We compute the gravitational magnification taking into account effects 
of the lensing and show that, whereas the innermost 7-ray spheres can be significantly 
magnified, there is little magnification either for very high 7-ray energies or for lower 
(radio) frequencies (because these wavelengths are emitted from larger regions). We 
analyze the temporal evolution of the gamma-ray magnification for sources moving in 
a caustic pattern field, where the combined effects of thousands of stars are taken into 
account using a numerical technique. We propose that some of the unidentified 7-ray 
sources (particularly some of those lying at high galactic latitude whose gamma-ray 
statistical properties are very similar to detected 7-ray blazars) are indeed the result of 
gravitational lensing magnification of background undetected Active Galactic Nuclei 
(AGNs). This is partly supported from a statistical point of view: we show herein as 
well, using the latest information from the Third EGRET Catalog, that high-latitude 
7-ray sources have similar averaged properties to already detected 7-ray AGNs. Some 
differences between both samples, regarding the mean flux level, could also be under- 
stood within the lensing model. With an adequate selection of lensing parameters, it 
is possible to explain a variety of 7-ray light curves with different time scales, includ- 
ing non- variable sources. The absence of strong radio counterparts could be naturally 
explained by differential magnification in the extended source formalism. The spectral 
evolution of the sources during microlensing events is calculated, revealing specific fea- 
tures that can be used to test the models with the next generation of both, orbital and 
ground-based 7-ray telescopes. 

Keywords: gamma-rays: observations, gamma-rays: theory, gravitational lensing, 
galaxies: active 

1 INTRODUCTION sky, remain unidentified at present. 



The first extragalactic 7-ray source detected was the quasar 
3C273, which was observed by the COS-B satellite in a par- 
ticularly active state in the f970s (Swanenburg et al. f978). 
Since then, many Active Galactic Nuclei (AGNs) have been 
detected at high energies, most of them belonging to the 
blazar class (e.g. Mukherjee 2001). The Third EGRET 
Catalog of point-like sources lists currently 67 detections 
labeled as AGNs (Hartman et al. 1999). Notwithstanding, 
a large number of 7-ray sources, scattered along the entire 



The unidentified 7-ray sources at low latitudes are 
probably related to several distinct galactic populations 
(Romero 2001). Among them there might be pulsars 
(e.g. Kaspi et al. 2000, Zhang et al. 2000, Torres et al. 
2001d; Camilo et al. 2001; DAmico et al. 2001), SNRs 
in interaction with molecular clouds (Combi & Romero 
1995, Sturner et al. 1996; Esposito et al. 1996; Combi et 
al. 1998; Combi et al. 2001, Butt et al. 2001, Torres et al. 
2002b), stellar-size black holes (Punsly 1998a, b; Punsly et 
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al. 2000), X-ray transients (Romero et al. 2001), persistent 
microquasars (Paredes et al. 2000, Kaufman Bernado et al. 
2002), and massive stars with strong stellar winds (Benaglia 
et al. 2001). Some of these kinds of stellar objects present 
statistical positional correlation with unidentified EGRET 
sources (far from what is expected as a random result, e.g. 
Romero et al. 1999, Torres et al. 2001b). Pulsars, however, 
remain as the only confirmed low-latitude population, since 
pulsed 7-ray radiation has been already detected for at 
least six different sources (Thompson et al. 1999, Thompson 
2001; Kaspi et al. 2000), five of them included in the Third 
EGRET Catalog (Hartman et al. 1999). 

Gehrels et al. (2000) have shown that the mid-latitude 
sources are different from the bright population of unidenti- 
fied sources along the Galactic plane. Some of the detections 
(5° < |6| < 30°) are thought to be associated with the 
Gould Belt (Grenier 2000, Gehrels et al. 2000), a starburst 
region lying at ~ 600 pc from Earth. Few other sources, 
at higher latitudes, could be the result of electrons being 
accelerated at the shock waves of forming clusters of 
galaxies (Totani & Kitayama 2000). However, for many of 
the high-latitude sources, no other explanation seems to 
be available than they are AGNs yet undetected at lower 
energies. This is particularly clear when one looks at the 
variability levels of the associated light curves: models 
requiring large acceleration region, like clusters of galaxies, 
would produce non-variable sources, contrary to what is 
found for most of the high latitude sources. 

All identified 67 EGRET AGNs are also strong radio 
sources with flat spectra, as expected from synchrotron 
jet-like sources where the 7-ray flux is the result of inverse 
Compton scattering (Mattox et al. 1997). However, no 
strong radio source appears within the contours of the 
unidentified high latitude EGRET sources. In this paper, 
we shall develop a model, briefly outlined by Torres et 
al. (2002a), which focuses precisely on that difference 
and provides an explanation of why some of the high 
latitude unidentified sources might not be detected at low 
frequencies. The main feature of such a model is that it 
will account for the 'j-ray properties of the high latitude 
gamma-ray detections resorting to differential gravitational 
lensing magnification of background, high-redshift, AGNs 
with otherwise undetected j-ray emission. Since these 
objects have different sizes at different wavelengths, differ- 
ential microlensing effects will lead to a magnification of 
the innermost ^-ray emitting regions, whereas the radio 
emission will be largely unmagnified, therefore remaining 
under the detection threshold. 

The gravitational light deflection effect by compact 
objects on background sources is commonly called mi- 
crolensing (e.g. Paczynski 1986). A source would be affected 
by different magnifications, depending on its position. 
Typically, source, lens, and observer move relative to each 
other, and therefore, this translates into a variable flux 
measured for the background source. Observationally, 
there are two interesting regimes of microlensing. Local 



microlensing deals with the light deflection effects by stars 
inside the Milky Way disk on stars in the Galactic bulge. 
Here the probability for a microlensing event is of order 
xlO" 6 . This means that it is necessary to monitor millions 
of stars in order to see a few occurrences. But despite this 
small probability, various teams have been very successful 
in detecting this kind of events in recent years (for a 
review, see Paczynski 1996). The other interesting regime 
of microlensing is usually called quasar microlensing, but it 
can be applied to any other compact source at moderate to 
high redshift. In this case, an intervening galaxy provides 
the surface mass density in stars (or other compact objects) 
which act as microlenses on the background quasar (for 
a review, see Wambsganss 2001). Recently, this kind of 
microlensing has been suggested for other astrophysical 
sources as well, e.g. gamma-ray bursts (Williams & Wijers 
1997), gamma-ray burst afterglows (Garnavich, Loeb & 
Stanek 2000; Mao & Loeb 2001; Koopmans & Wambsganss 
2001), and superluminal shocks in extragalactic radio 
sources (Romero et al 1995, Koopmans & de Bruyn 2000). 
The gamma-ray sources discussed here are another type of 
astrophysical objects for which microlensing possibly plays 
an important role. 

We recall that gravitational light deflection is basically 
an achromatic phenomenon, being a geometric effect pre- 
dicted by General Relativity, i.e. the deflection angle does 
not depend on the energy of the photon. However, it is 
nevertheless possible to have chromaticity effects when the 
size of the source changes with the observing wavelength. A 
large source is typically less affected by a microlensing mag- 
nification than a small source (Wambsganss & Paczynski 
1991). This size-induced chromaticity will be an essential 
ingredient of our model. 

The structure of the paper is as follows. In Section 2 we 
present a statistical study of the main characteristics of the 
sample of unidentified 7-ray sources at high latitudes, in- 
cluding variability. In Section 3 we review the formalism for 
extragalactic gravitational microlensing in both the point- 
like and extended source cases. Section 4 comments on the 
inner structure of 7-ray emitting AGNs considered here as 
background sources. Section 5 contains our results for the 
flux and spectral evolution of single microlensing events in a 
variety of situations and different types of host galaxies. Sec- 
tion 6 deals with the optical depth problem and the expected 
number of microlensing events. Section 7 presents our results 
for full numerical modeling of caustic patterns for different 
galactic lensing parameters. Precise predictions for the light 
curves at different frequencies are also shown there. We fi- 
nally close with some concluding remarks in Section 8. 



2 UNIDENTIFIED 7-RAY SOURCES AT HIGH 
LATITUDES 

2.1 Sample and photon spectral index 

Previous population studies using the Second EGRET 
Catalog have already remarked that part of the sample 
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of high latitude unidentified sources is consistent with an 
isotropic population, a fact that supports an extragalactic 
origin for these detections (Ozel & Thompson 1996). In 
what follows, we shall make a comparison between the 
properties of identified 7-ray AGNs in the Third EGRET 
Catalog and high-latitude unidentified sources. We shall 
choose the lower cut-off in latitude as |6| = 30°, in order 
to avoid possible contamination from Gould Belt sources. 
There are 45 3EG unidentified sources within this latitude 
range; we provide details on these sources in Table 

Figure ^ shows the distribution of the 7-ray photon 
spectral index for both sets of sources, 45 unidentified 
EGRET sources and 67 detected EGRET AGNs. The mean 
value of the photon index is 2.36±0.36 for AGNs, and 
2.49±0.34 for the unidentified detections. They are compat- 
ible within the uncertainties and, on average, steeper than 
what is observed for low latitude sources, which are thought 
to belong to our own Galaxy. 
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2.2 Variability 

The study of the time variability of 7-ray sources is of 
fundamental importance. Several models for 7-ray sources 
in our Galaxy predict non- variable emission during the time 
scale of EGRET observations. AGNs, on the contrary, are 
expected to present a variable flux emission. Variability 
is, in case we were able to quantify it with some degree 
of precision, a powerful tool to probe the nature of the 
sources. Visual inspection of the flux evolution through 
the different viewing periods is obviously a first indication 
of the variability status of any given source. However, 
fluxes are usually the result of only a handful of incoming 
photons, experimental errors are sometimes huge and their 
origin uncertain, and consequently more reliable ways of 
quantizing the flux evolution should be devised: these are 
known as variability indices. Two such indices have been 
recently introduced in the literature and applied to 3EG 
sources so far (Tompkins 1999 (index r), Torres et al. 2001a 
(index /)). In general, statistical results from these two 
indices are well correlated (see Torres et al. 2001c for a 
discussion). Here we shall adopt the index /, previously 
used in blazar variability analysis (Romero et al. 1995) 
and applied to some of the 3EG sources by Zhang et al. 
(2000) and Torres et al. (2001a,c) as our main quantitative 
evaluation of variability, although the results for r are also 
given in Table 1. The basic idea behind the index / is to 
directly compare the flux variation of any given source with 
that shown by known 7-ray pulsars, which are assumed 
to be an intrinsically non-variable population. This index, 
contrary to Tompkins' index r (Tompkins 1999), uses only 
the publicly available data of the 3EG Catalog. 

Let us recall the basic elements that are used to define 
the /-index. Firstly, a mean weighted value for the EGRET 
flux is computed: 



where N vp is the number of single viewing periods, F(i) the 
observed flux in the i th -period, and e(i) the corresponding 
error. For those observations in which the significance (yTS 
in the EGRET catalog) is greater than 2a, the error is e(i) = 
F(i) /\/TS. For those observations which are in fact upper 
bounds on the flux, it is assumed that both F(i) and e(i) are 
half the value of the upper bound. The fluctuation index fj, 
is defined as: 



/i = 100 x a sd x (F) 1 



(2) 



In this expression, <r s d is the standard deviation of the flux 
measurements. This fluctuation index is also computed for 
the confirmed 7-ray pulsars in the 3EG catalog, assuming the 
physical criterion that pulsars are non-variable 7-ray sources. 
The averaged statistical index of variability, I, is then given 
by the ratio 

j ^source ^\ 

< /I > pulsars 

Once the index is defined, we need to clarify the thresholds 
for variability. Following Torres et al. (2001c), clearly 
variable sources will be those with I > 5, possibly variable 
sources will have 2.5 < / < 5, non- variable sources will 
have I < 1.7, and the remaining sources will be considered 
as dubious cases. These are very conservative cut-offs: 
I > 5 means that we are asking for the value of I to be 8a 
above that of pulsars in order to classify a source as variable. 

In Figure ^ we compare the /-index distribution for 
the samples under analysis. The mean value for AGNs (left 
panel) is 3.3±2.6. A possible peak in the plot is seen at 
/ = 2.5, which represents a value 4a above that presented 
by pulsars. An apparently extra peak appears at / ~ 6. 
Clearly, most of the AGNs are likely variable sources. The 
mean for the unidentified sources (Figure ^, right panel) is 
also high: 3.0±2.3. There are, again, apparently two distinct 
peaks in the distribution, located at / ~ 2 and / ~ 6. 

In Figure ^ (left panel) we show the variability index / 
versus the Galactic latitude. The constraint we are imposing 
on our sample of unidentified sources (to have |6| > 30°) can 
be clearly noticed in the bottom plot. There is not a clear 
dependence of the variability index with latitude, neither for 
AGNs nor for unidentified sources. The same happens in the 
plots of Figure ^| (right panel) , where we show the variability 
index versus the photon spectral index. An apparent trend 
of increasing the variability status for the steepest sources, 
already noticed by Torres et al. (2001a) and Reimer (2001), 
is shown in this figure. However, this is not conclusive since 
results for a Spearman Rank test are in the range of a few 
percent for this to be a random phenomenon. An overall 
characteristic of Figures ^| and ^ is that both samples look 
quite similar, with no apparent strong deviation from each 
other shown in terms of variability or photon spectral index 
distributions. 
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Table 1. The 45 unidentified sources considered in the analysis. We list their 3EG Catalog name, their Galactic coordinates, spectral 
index, variability index, and the values of (F) used to define / in Eq. (g). We also provide the 3EG P1234 fluxes, F. (F) and F are in units 
of 10 — 8 photons cm~ 2 s _1 . The columns labeled t, T m i„, and r max , give the central value of Tompkins' (1999) index for variability and 
their 68% CL lower and upper limit deviations, respectively. The pulsar population has (t) < 0.1, whereas typical AGNs have (r) ~ 0.7. 
Extreme upper limits for r, whose maximum is 10000, imply possible strong variability. 
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0245+1758 


157.62 


-37 


.11 


2.61 


2.74 


2.63 


0.73 


2287 


11.3 


8.8 


0404+0700 


184.00 


-32 


.15 


2.65 


1.50 


0.34 


0.00 


1.65 


13.5 


11.1 


0512-6150 


271.25 


-35 


.28 


2.40 


2.34 


0.00 


0.00 


0.55 


10.8 


7.2 


0530-3626 


240.94 


-31 


.29 


2.63 


1.62 


17.8 


15.8 


0.61 


0.15 


2.28 


0808+4844 


170.46 


32 


.48 


2.15 


0.90 


0.00 


0.00 


0.39 


12.2 


10.7 


0808+5114 


167.51 


32 


.66 


2.76 


1.53 


0.00 


0.00 


0.73 


13.5 


8.7 


0910+6556 


148.30 


38 


.56 


2.20 


1.79 


0.49 


0.00 


1.14 


9.2 


5.9 


1457-1903 


339.88 


34 


.60 


2.67 


2.72 


0.42 


0.00 


3.64 


13.9 


8.1 


1504-1537 


344.04 


36 


.38 


- 


2.73 


10.33 


1.21 


9999. 


12.9 


8.8 


1600-0351 


6.30 


34 


.81 


2.65 


7.26 


68.17 


0.00 


9999. 


9.7 


9.9 


1621+8203 


115.53 


31 


.77 


2.29 


0.34 


0.00 


0.00 


0.29 


11.5 


7.4 


1733+6017 


89.12 


32 


.94 


3.00 


1.82 


0.39 


0.00 


1.38 


16.9 


8.7 


1958-4443 


354.85 


-30 


.13 


- 


7.43 


58.02 


5.85 


9999. 


11.3 


6.4 


2034-3110 


12.25 


-34 


.64 


3.43 


5.26 


2.88 


0.89 


155. 


6.8 


5.2 


2219-7941 


310.64 


-35 


.06 


2.50 


1.37 


0.00 


0.00 


0.51 


19.8 


13.5 


2243+1509 


82.69 


-37 


.49 


- 


11.11 


3.42 


0.88 


3097 


7.4 


9.9 


2248+1745 


86.00 


-36 


.17 


2.11 


2.20 


1.07 


0.43 


3.98 


20.2 


12.9 


2255+1943 


89.03 


-35 


.43 


2.36 


5.54 


2.31 


0.80 


48.6 


14.2 


5.8 


0038-0949 


112.69 


-72 


.44 


2.70 


2.71 


0.00 


0.00 


0.89 


15.3 


12.0 


0118+0248 


136.23 


-59 


.36 


2.63 


2.28 


5.17 


0.90 


9999. 


11.5 


5.1 


0130-1758 


169.71 


-77 


.11 


2.50 


0.10 


0.00 


0.00 


0.38 


13.5 


11.6 


0159-3603 


248.89 


-73 


.04 


2.89 


1.26 


0.00 


0.00 


1.16 


11.9 


9.8 


0215+1123 


153.75 


-46 


.37 


2.03 


3.67 


10.06 


1.19 


9999. 


8.0 


9.3 


0253-0345 


179.70 


-52 


.56 


- 


7.99 


16.44 


1.38 


9999. 


5.0 


6.2 


0348-5708 


269.35 


-46 


.79 


- 


4.76 


6.60 


1.29 


9999. 


5.8 


3.8 


0917+4427 


176.11 


14 


.19 


2.19 


1.69 


0.00 


0.00 


0.34 


18.6 


13.8 


1009+4855 


166.87 


51 


.99 


1.90 


1.14 


0.00 


0.00 


0.60 


7.7 


4.8 


1052+5718 


149.47 


53 


.27 


2.51 


2.02 


0.21 


0.00 


0.74 


6.8 


5.0 


1133+0033 


264.52 


57 


.48 


2.73 


4.44 


0.71 


0.16 


2.00 


9.1 


3.7 


1134-1530 


277.04 


43 


.48 


2.70 


3.35 


2.85 


1.11 


51.5 


17.2 


9.9 


1212+2304 


235.57 


80 


.32 


2.76 


5.65 


78.82 


0.00 


9999. 


6.3 


3.3 


1219-1520 


291.56 


46 


.82 


2.52 


3.16 


1.78 


0.74 


13.7 


8.9 


4.1 


1222+2315 


241.87 


82 


.39 




1.95 


61.09 


2.44 


9999. 


6.9 


5.7 


1227+4302 


138.63 


73 


.33 




3.20 


61.09 


2.44 


9999. 


6.7 


4.6 


1234-1318 


296.43 


49 


.34 


2.09 


2.03 


0.42 


0.12 


0.81 


11.0 


7.3 


1235+0233 


293.28 


65 


.13 


2.39 


1.04 


0.23 


0.00 


0.65 


10.5 


6.8 


1236+0457 


292.59 


67 


.52 


2.48 


1.70 


0.00 


0.00 


1.45 


7.3 


6.5 


1310-0517 


311.69 


57 


.25 


2.34 


1.28 


2.94 


1.69 


7.92 


11.4 


7.9 


1323+2200 


359.33 


81 


.15 


1.86 


5.17 


2.69 


0.93 


46.8 


10.1 


5.2 


1337+5029 


105.40 


65 


.01 


1.83 


2.85 


0.54 


0.00 


1.35 


10.0 


9.2 


1347+2932 


47.31 


77 


.50 


2.51 


1.10 


0.48 


0.00 


1.45 


15.3 


9.6 


1424+3734 


66.82 


67 


.76 


3.25 


1.90 


0.01 


0.00 


9999. 


18.0 




2241-6736 


319.81 


-45 


.02 


2.39 


1.25 


0.00 


0.00 


1.09 


16.6 




2251-1341 


52.48 


-58 


.91 


2.43 


5.17 


9.49 


1.58 


9999. 


9.7 


6.5 


2255-5012 


338.75 


-58 


.12 


2.79 


1.59 


0.41 


0.00 


1.46 


12.7 


9.2 



2.3 Fluxes and possible radio counterparts 

In Figure ^ we show the EGRET averaged flux as a function 
of the photon spectral index and the variability index, 
respectively. We notice that, although there is no apparent 
difference in the form of the distribution for both samples 
under consideration, there is a clear contrast on the flux 
values: Whereas most identified 7-ray AGNs have fluxes 



above 10 -7 photons cm -2 s _1 , most of the unidentified 
sources present lower values. This is consistent with what 
was presented by Gehrels et al. (2000) for sources at 
latitudes |6| > 5°. 

This difference in the flux values is also translated into 
the Log N - Log F plots we present in Figure pi It can 
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Figure 1. Photon spectral index comparison. The left panel shows the distribution for the 67 detected AGNs, dubbed A, in the Third 
EGRET Catalog. The right panel shows the corresponding distribution for the |b| > 30° unidentified sources. 



be seen there that the linear fits differ significantly. AGNs 
present a fit close to what is expected for an isotropic 
and uniform population (F~ 3 ^ 2 ), and also similar to what 
was found using the 2EG Catalog (Ozel & Thompson 
1996). The unidentified sources, however, present a steeper 
dependence. The difference in the flux levels is also shown 
in the a;-axis. However, the analysis of the result for the 
sample of unidentified sources should be done with extra 
care, since the errors are far larger, as well as the number of 
sources considered is smaller. Additionally, AGNs present 
an apparent lack of sources at F ~ 30 x 1CP 8 photons 
cm -2 s _1 , which should be confirmed or falsified by future 
observations. 

Reimer & Thompson (2001) studied in detail the log N 
- log F plots obtained from 3EG sources, but including also 
those sources with lower confidence level (which did not 
appear in the published version of the 3EG catalog). They 
found that there is a very pronounced contrast between 
average and peak flux representation in a log N - log F 
diagram for the sources above |fc| > 30°. This is due to the 
fact that sources at high latitudes are mostly detected only 
in some (or in many cases, only in one) viewing periods 
(see below), when they show their peak flux, leaving the 
average over the four phases of the experiment in a much 
lower value. The differences in fluxes between the peak 



detections of both distributions, although still present, are 
not so strong as the ones presented in the P1234 averaged 
values. As we shall see below, sources showing large fluxes 
only in one viewing period could be particularly suitable to 
be explained by microlensing of gamma-ray blazars. 

The differences in the Log N - Log F plots can be point- 
ing towards one of the two following possibilities: 

(i) We are looking at (at least) two different populations; 
for instance, AGNs and a new halo class of high-energy ob- 
jects. 

(ii) These samples are formed mostly by the same kind 
of objects (AGNs) but they present different 7-ray flux lev- 
els. This difference could be produced as an extrinsic effect 
when the sources are farther away from us than those which 
produce the most energetic detections. 

In this paper we shall explore a plausible situation in 
which the second possibility prevails. If behind both samples 
there is actually a single class of extragalactic objects, the 
combined Log N - Log F plot should approximately follow 
a N oc F~ 3 / 2 law. The result, given in Figure ^| confirms 
this and it is close to the isotropic and uniform expectation 
(see also the comments by Reimer & Thompson 2001). 
This seems to be suggesting that the sample of unidentified 
sources under consideration is formed by many weak AGNs, 



© 0000 RAS, MNRAS 000, 000-000 



6 Torres, Romero, Eiroa, Wambsganss, & Pessah 




Figure 2. Variability index comparison. The left panel shows the distribution for the 67 detected AGNs, dubbed A, in the Third EGRET 
Catalog. The right panel shows the corresponding distribution for the |f>| > 30° unidentified sources. 



which at the same time presents a low (i.e. under the 
thresholds of the corresponding surveys) radio emission. We 
note that the apparent bump at F ~ 30 x 10~ 8 photons 
cm -2 s _1 continues to appear in the combined plot, since no 
unidentified sources present such relatively high 7-ray fluxes. 

Recently, Mattox et al. (2001) have presented a quanti- 
tative analysis of potential radio identifications for all 3EG 
sources. They used radio surveys at 5 GHz, as it was done 
previously for the 2EG Catalog (Mattox et al. 1997), and 
evaluated an a priori probability for these associations to be 
physical, based on the positional offsets and radio fluxes of 
the proposed counterparts. They found that 45 out of the 67 
3EG sources classified with 'A' by Hartman et al. (1999) were 
among the EGRET identifications with the highest proba- 
bility of being correct. Only one extra possible association in 
the list of these most likely identifications was not dubbed 
'A' in the 3EG Catalog. For each of these 46 associations, 
they have compiled radio fluxes at 5 GHz, and when avail- 
able, also those obtained with VLBI. In Figure ^ we show the 
7-ray flux of each of these 3EG sources (note that we plot 
the P1234 EGRET flux of the 3EG source from Hartman et 
al. 1999, not the flux of the AGN quoted, for instance, in 
Mattox et al. 1997) as a function of the radio flux at 5 GHz 
of the likely counterparts. Most of the detections present a 
radio flux above 1 Jy, and there is an apparent trend to be- 



come more radio loud when the observed 7-ray flux is higher. 
This is the expected behaviour when the emitted 7-rays have 
their origin in inverse Compton interactions of the same par- 
ticle population (leptons) that generates the radio emission, 
targeting a soft photon field. We have superposed a linear 
function to the radio-gamma data that is also shown in Fig- 
ure There, we indicate with a dashed line, the extension 
of this linear fit to the region where there are no 3EG 'A' 
AGN sources. The big vertical box signals the threshold for 
detectability in the radio surveys used to search for counter- 
parts (~ 30 mjy). The middle, lighter colored box, signals 
the range of 7-ray fluxes for the unidentified sources consid- 
ered in this paper. It is apparent, then, that if weak AGNs 
were to approximately follow the linear fit, we could find 
several 7-ray sources without significant radio flux. Many of 
them could be those unidentified sources we are studying 
here. In addition, if 7-ray sources are affected by a differen- 
tial gravitational lensing effect, this process, as we shall show 
below, would enhance only the 7-ray emission, keeping the 
radio fluxes at low levels. This mechanism, then, would be 
in agreement with what is shown in Figure [?], provided the 
associated sources are within the middle box on the left. Of 
course, this cannot apply to all unidentified sources because, 
otherwise, it would result in a hole in the source distribu- 
tion between the already detected AGNs and the candidates, 
at radio flux levels of ~ 100 mjy. It should be remembered, 
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Figure 3. Variability index / (Torres et al. 2001) versus Galactic latitude (left) and versus photon spectral index (right). 



anyhow, that sources at high latitude are preferentially iden- 
tified by their peak flux, which can be much higher than the 
average. Although population studies are a powerful tool to 
study the nature of the unidentified detections, they should 
be supplemented by a source-by-source analysis. 



2.4 Light curves 

The most interesting candidates in order to pursue further 
study seem to be those sources presenting the highest lev- 
els of variability. Most AGNs are, as we have already said, 
variable sources, even with 7-ray emission fluctuating from 
detection to non-detection in successive viewing periods. We 
show some examples of typical light curves in Figure ^. This 
kind of light curves are also shown by many unidentified 
sources, what is consistent with the high levels of variability 
they present. 



2.5 Assessment 

To finish this section we remark our main conclusion up to 
here, namely, that it is likely that some of the high-latitude 
unidentified sources be not more than otherwise undetected 
AGNs, presenting a low or nil (below any current detection 
threshold) radio flux. General gamma-ray characteristics of 
both 7-ray blazars and high-latitude unidentified sources are 
very similar. There remains, however, the question of why, 
whereas most 7-loud blazars present radio flux at the Jy 
level, the unidentified sources have no strong radio coun- 
terpart at all. If the same mechanism for 7-ray production 
operates in both groups of objects, why are they so different 
at lower energies? We shall argue in the next sections that 
extrinsic effects can result in such a behaviour. 



3 MICROLENSING 

3.1 Point-like source and point-like lens 

Let us consider a background and weak 7-ray emitting 
blazar, whose GeV flux is well below the EGRET sensitiv- 
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Figure 4. Source fluxes (P1234, Hartman et al. 1999) versus photon spectral index (left) and variability index (right). 



ity threshold and whose radio flux is at the mjy level. We 
shall first use the Chang & Refsdal (1984) model. Assuming 
a galaxy is interposed in the line of sight, the lens equation 
in the lens plane for a point source is (e.g. Schneider et al. 
1992): 



r - r - R E 



K + 7 
K 







r - do = 0, 



(4) 



where the coordinate system is centered on the microlens (a 
star within the galaxy) , with the orientation of the orthonor- 
mal basis {ei,e2} chosen to diagonalize the quadrupole ma- 
trix; the source is at ro and the image position is r. The 
third term in equation (^) arises from the deflection in the 
lens plane due to the microlens, here considered as a point 
mass M. Re is the usual Einstein radius 



Re 



AGM DoiD h 



(5) 



with D os the observer-source distance, D a i the observer-lens 
distance and Di s the lens-source distance. The fourth and 
fifth terms in Eq. M) arise from the deflection imprinted by 



the galaxy as a whole. In the fourth term, k and 7 are the 
focusing and the shear of the galaxy, at the lens position, 
respectively, do depends on the deflection imprinted by the 
galaxy as a whole at the location of the microlens; its only 
effect is to change the unperturbed source position ro by a 
constant. We shall ignore do, assuming a displaced source 
position s = ro + do in the lens plane. 

Defining new coordinates X and Y, in the lens and the 
source plane respectively, as 



X = 
and 
Y = 



y/\l-K + y\ 
Re ' 



Re\/\1 - k + 7| 
the lens equation becomes 



1 



IXI' 



(6) 
(7) 

(8) 



where 
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Figure 5. Log N-Log F comparison. The value of the linear fit slope (and error within parentheses) is shown for each case. F is the 3EG 
P1234 flux. 



s = sign (1 — n + 7) , (9) 
and 

A=^-X (10) 

1 — K + 7 

The solution of equation (^) can be found reducing the prob- 
lem to a fourth order equation for X 2 (Schneider et al. 1992). 
In order to do it we introduce new coordinates 



X\ = X cos <j>, 
X2 = X sin cj>, 
so we obtain 

Y 1 = eAX 1 -^ = {eA- ±)xi 



and 



X 2 



Y 2 = eX 2 - = e - — )X sin 



1 



X 2 V X 2 , 
Solving for cos 4> and for sin <f>, we get 
Yi 



and 



X(eA- 1/X 2 )' 



Y 2 



r X(e-l/X 2 )' 
© 0000 RAS, MNRAS 000, 000-000 



(11) 
(12) 

(13) 

(14) 

(15) 

(16) 



This readily implies 



= 1, 



X 2 (eA-l/X 2 ) 2 ' X 2 (e-1/X 2 ) 2 
that translates into 

A 2 X s - [2eA (A + 1) + if + A 2 Y 2 2 ] X 6 + 
[A 2 + 4A + 1 + 2e (Y 2 + AY 2 2 )] X 4 - 
[2e (A + 1) + Y'l + Yi] X 2 + 1 = 0. 



(17) 



(18) 



Here, X = |X| and the sub-indices stand for the different 
components of the respective vectors. Equation (^) can be 
solved by any standard method. Since it is a fourth order 
equation, we can have zero, two or four real solutions to the 
lens equation. The number of images is then given by the 
number of real and positive solutions of this equation. For 
any real and positive solution, X 2 , we can take the positive 
value of X and obtain the position of the image using the 
definitions given in Eqs. (|l5|-|l^) and (p]-[l^) above. 

The magnification A = I b s /Io for any image is in turn 
given by the determinant of the Jacobian matrix. This can 
be computed using the chain rule for derivatives 

ds_dY_ _ ds dY dX 
dYlfr ~ dYdX^lh' 
that finally yields 



(19) 
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Figure 8. Examples of 7-ray light curves of 3EG 'A' AGNs (upper panel) and unidentified sources with |6| > 30° (lower panel). Cases 
chosen among the most variable ones are shown here. Points with arrows represent upper limits on the flux. Spline curves do not represent 



actual fits. 
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Figure 6. Log N-Log F for the combined samples. The value of 
the linear fit slope (and error within parentheses) is shown. 



A = 

Since 

dY 

9X = 

then 

A = 



det 



(sr) 



1-/S + 7 



eA 



-r — jjr- 

2x1X2 



\1-K- 



2X1X2 
X4 

x^ - 



(20) 



(21) 



I A (X? + X*) A + g (1 — A) (X? - Xf ) - 1| 
The total magnification for a point source is 



A = ^ Ai 



(22) 



(23) 



where n is the number of images and Ai are obtained 
replacing the solutions of Eq. (^) in Eq. (p2|). 

The curves for which the determinant of the Jacobian 
matrix is zero are called critical curves. For the Chang- 
Resfdal lens, these critical curves are Cassini ovals given by 
the equation 



A (Xl + Xl) + e (1 - A) {x\ - X\) -1 = 0. 



(24) 



Figure 7. Radio and 7-ray fluxes of the 46 most likely AGN 
detections. The big vertical box signals the threshold sensibility in 
the radio surveys used to search for counterparts (~ 30 mjy). The 
middle box signals the range of 7-ray fluxes for the unidentified 
sources considered in these paper. 



The caustics are the mappings of the critical curves, through 
the lens equation, onto the source plane. When the source 
crosses a caustic, the number of images changes by two, 
and, if a point source is on a caustic, the magnification 
diverges. In realistic situations, sources are extended instead 
of point-like, and when the source is on (or near) a caustic 
the magnification is finite (although large). 

When 7 = (no shear case), A = 1, and the lens geom- 
etry becomes axially symmetric (notice also that the caus- 
tic curve collapses to a point). Then, we can use the one- 
dimensional lens equation: 



Y = eX- 



1 



X ' 

which has as solutions 



x± = - ± Vx +£ 



(25) 



(26) 



for Y 2 > — 4e. Hence, for e = 1 we always have two solutions, 
and when £ = — 1 there are two solutions for \Y\ > 2 and 
none for \Y\ < 2. The magnification of any image is 
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A ± = 



l 



X 



± 



and the total magnification is given by 
A = A+ + A-. 



(27) 



(28) 

Replacing Eq. ( [27] ) in Eq. ( |28[ ) and using Eq. (Q, we arrive 
at 

V 2 + 2e 



A 



1 - k\ \Y\ VY 2 + 4e ' 



(29) 



Some tricky -although straightforward- algebra used to ar- 
rive at the previous equation is described in the Appendix. 
In this case, we have the critical curve at X = 1, which 
gives a degenerate caustic at Y — for e = 1, and a circular 
caustic of radius Y = 2 for e = — 1. 



3.2 Extended source 

For an extended circular source, the light curve is given by 
(see e.g. Han et al. 2000): 

(■27T ff. 



So * Jo 3 J ( r ' ( r ' <fi) rdrd V 

ST Jo" J ( r ' t p) rdrdL P 



(30) 



where (r, y>) are polar coordinates in a reference frame placed 
in the center of the source, r s is the radius of the source and 
T(r, cp) is the surface intensity distribution of the source. For 
a uniform distribution this produces a simplified expression, 

f2w nr. 



A 



Jo* Jo' A {r,<p)rdrd<p 



(31) 



Introducing the dimensionless radial coordinate R = r/Rm, 
we obtain 



S? Io Rs A (R,p)RdRd<p 



(32) 



where R s = r s /RE is the corresponding dimensionless 
radius of the source. 

We assume that the microlens is moving with constant 
velocity v and we choose the origin of time (t — 0) as the 
instant of closest approach between the lens and the source. 
Then, if the center of the source is placed at t = in b = 
(61,62), the position of any point of the source with polar 
coordinates (r, ip) is: 



si(t) = 61 — vt cos 9 + r cos ip, 
S2 (t) — 62 — vt sin 8 + r sin ip, 



(33) 
(34) 



where v = jvj , 6 is the angle between v and ei, < r < r s 
and < ip < 2-7T. Figure [] shows these geometrical consider- 
ations. 



Yl = 
Y 2 = 



Plugging Eqs. into Eq. (§), we obtain 

61 — vt cos 8 + r cos ip 
Re\J\1 - « + 7| 

62 — vt sin 8 + r sin ip 



(35) 
(36) 



Source 




Lens 



Figure 9. Geometrical sketch of an extended source and physical 
variables. All points of the source must be taken into account 
when computing the magnification in a microlensing event. 



In units of the Einstein radius, these equations transform 
into 

Bi — Tcosd + Rcosip 



Yi 



Y 2 = 



B 2 - T sin 8 + R sin ip 

V|l-« + 7l 



(37) 



(38) 



where T = vt/R E , Bi = bi/R B and B 2 = b 2 /R E ■ When 
7 = (no shear), we can take 8 = and B = (0, Bo) in the 
above expressions without loosing generality. 



4 7- RAY BLAZARS AS SOURCES 

The fact that some 7-ray blazars have been observed to 
flare dramatically on time scales of days imposes severe con- 
straints on the size of the emitting region. The optical depth 
for intrinsic 7 + 7 — > e + +e~ attenuation is (e.g. Schlickeiser 
1996): 

(IT 



-I. 



(39) 



47rc (e) 

where or is the Thompson cross section, n 7 the 7-ray photon 
density, R < ct v the source size inferred from the intrinsic 
variability time scale, < e > the mean photon energy, and 
I the compactness parameter defined as the ratio of the in- 
trinsic source luminosity to its radius. The optical depth can 
be written as: 
£48 



r > 200 



t v /l day ' 



(40) 
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with 1/48 the luminosity in units of 10 48 erg s _1 . For typical 
values t y ~ 1 day and L48 ~ 1, the source is opaque, 
contrary to the observed fact that 7-ray blazars present a 
power-law 7-ray spectrum over several decades of energy. 
This rules out isotropic emission in the rest frame. The 
emission, consequently, should be beamed. It is usually 
thought to be produced in a relativistic jet through inverse 
Compton scattering of lower energy photons (e.g. Blandford 
& Levinson 1995). The soft, seed photons for the Inverse 
Compton process could originate as synchrotron emission 
from within the jet, or they could come from the accretion 
disk surrounding the central supermassive compact object, 
or they could be disk radiation reprocessed in the broad line 
region. In the last two cases the seed photons are external 
to the jet itself. 

In addition to these leptonic models, some hadronic 
alternatives have been proposed in the literature. The 7-ray 
emission would be produced in this case by relativistic 
protons interacting with ambient matter, radiation fields, or 
the magnetic field of the jet. For reviews and references the 
reader can see von Montigny et al. (1995) and Mukherjee 
(2001). 

Independently of how the 7-rays are produced, they 
must traverse the strong X-ray field produced in the 
innermost region of the accretion disk. The observed 7-ray 
photons cannot originate from a too small radius, otherwise 
they will be absorbed through pair creation in the disk 
photosphere (e.g. Becker & Kafatos 1995, Blandford & 
Levinson 1995). This naturally leads to the concept of 
7-spheres in AGNs: for each 7-ray photon energy there is 
a radius r 7 beyond which the pair production opacity to 
infinity equals unity (Blandford & Lenvinson 1995). The 
size of the 7-sphere will depend on both the energy of the 
7-ray photons and the soft photon flux. 



where R~,{E) — t 1 (E)/Re and R re t = r Te {/RE- 

We shall assume that the intensity of the source (with- 
out yet being lensed) is uniform, and that its spectrum ap- 
proximately follows a power law 

Io(E) oc E~ ( , 

with £ G (1.7,2.7) (Krolik, 1999). Then 

-t 



(43) 



Io(E) = I„ 



( E le f ) 



(44) 



where J rG f is the intensity of the reference source. The surface 
intensity distribution of the source finally ends being 



ME) = 



HE) 



n(R,(E)) 2 - 

Since To {E) does not depend on (R, <p) we can write 



A = 



Jo Jo Aq(R, tp)RdRdtp 



niR^E)) 2 
Using that A = I/Io and Eq. (Q, we have 

/ = AIo = AI Ie{ ( rf-) \ 

We define, for plotting purposes, the ratio 

/ . / E n -e 



Itel ( E TO { 



(45) 



(46) 



(47) 



(48) 



as the intensity in units of I Ic f. We shall as well define 
the dimensionless impact parameter u — B/R^, with B — 
■J Bf + B\. Throughout the rest of this paper we shall adopt 
p = 1.1 , £ = 2, and a reference source with dimensionless 
radius R rc { = 2 x 10~ 3 at a 7-ray energy of i? re f = O.lGeV. 
Results are, however, of a generic nature, and remain self- 
similar under reasonable variations of these parameters. 



For an isotropic, power-law, central source of soft pho- 
tons scattered by free electrons in a warped disk, Blandford 
& Levinson (1995) obtain: 



r 7 (£) tx E p , 



(41) 



with p depending on the details of the central source. A 
similar result is obtained for pure disk emission (Becker & 
Kafatos 1995, Romero et al. 2000). Typically, p £ [1, 2]. The 
larger 7-spheres, then, are those for the higher photon ener- 
gies. This energy-dependency of the source size will naturally 
lead to chromaticity effects during microlensing events. 



5 LIGHT CURVES AND SPECTRA OF 

SINGLE POINT LENSES EMBEDDED IN 
SMOOTHLY DISTRIBUTED MATTER 

Now we focus our study on the lensing of an extragalactic 
7-ray source. For immediate use, we shall define a reference 
source having a radius r re f and a 7-ray energy E le f such that 



R~z (E) — R FI 



V E Te f ) 



(42) 



5.1 Galactic model 

We shall now assume a parameterization for the mass distri- 
bution in the interposed galaxy. This will allow us to com- 
pute the values of k s and 7. In this Section, we assume for 
illustrative purposes that all matter in the lensing galaxy 
is smoothly distributed except for a single point lens. The 
full non-linear treatment for large number of lenses will be 
done in Section 7. If the center of the galaxy is located at 
x c = (b, 0) in the lens plane, and if its total mass is described 
by a surface density E given by 



(49) 



where E c is the density at the center of the galaxy, a g is 
the core radius of the galaxy, and d is the distance of the 
microlens to the center of the galaxy, then the focusing and 
shear can be written as (Schneider et al. 1992, Romero et al. 
1995): 



S Cl-it y/(lg +d 2 ' 



(50) 
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„_ S e / a g 2a 3 ^al + d? 2o|\ 

where E cr it is the critical density, defined by 



c 2 D s 



A-kGDiDu 

Defining m = E c /E cr it and \ — d/a g , we have 
1 



7 = m 



V 1 + X 2 

( 1 _ 2y^ + ^ 




When x "C 1, 7 adopts a simpler expression, 

.2 ,3 



7 = m 



(52) 



(53) 



(54) 



(55) 



If the microlens is placed at a distance from the center of 
the galaxy far smaller than the core radius, d <C a g , then 
X ~ and the lensing parameters result in 7 = 0, k 3 = m. 
We shall then adopt 7 = 0, k s = 0.8, in the zero shear case, 
as an example. To study other situations, where 7 can not 
be neglected, we shall take, m = 2, \ = 10 =>■ k s = 0.20, 
7 = -0.16.. 

5.2 Model results and analysis for a single point 
lens 

Results for the first group of parameters (shear equal to 
zero) are plotted in Figure h], for different impact distances. 
It can be seen there that background AGN "/-ray emission 
above 100 MeV can be magnified by a very large factor in 
single events. This magnification can make an otherwise 
unnoticed background source to exceed the detection 
threshold, and to be detected only when lensed. Note also 
the chromaticity effect: the 7 -sphere corresponding to 10 
GeV, whose size is similar to the radio emitting regions 
of AGNs, is negligibly magnified, while the lower energy 
curves -but still with energies in the EGRET range, 100-500 
MeV- all show significant magnifications This phenomenon 
also has a spectral signature, produced by the differential 
magnification of the different 7-spheres (see the discussion 
below). The magnification grows with the focusing. The 
same happens, as expected, when the impact parameter is 
smaller. 

A critical requirement for such a microlensing event to 
occur is that the size of the background source projected 
onto the lens plane is not larger than the Einstein ring of 
the lensing mass. Only background sources whose size is 
a fraction of the Einstein radius will then be significantly 
magnified. Since AGNs have emission regions of different 
scales for different radiation wavelengths, we expect a 
differential magnification, as observed in Figure llfj. The 
most internal regions of AGNs, responsible for the 7-ray 
emission, have sizes about 10 14 -10 15 cm (Blandford & 
Levinson 1995). Then, they can be significantly magnified. 
At the same time, since radio emission is originated far 



down the jet, the sizes of its emitting region (> 10 17 cm) 
exceeds the Einstein radius of the lenses and lead to the 
absence of counterparts at lower frequencies. 

In Figure [hJ we present a completely different behavior 
produced when the lensing parameter 7, the shear, is 
different from zero. We show the results for n 3 = 0.20, 
7 = —0.16, and 6 — 0°. 6 (see Figure is just the angle 
between the direction of the lens movement and the ei unit 
vector, in the source plane. 

In the cases of Figure [n], for the less energetic 7-spheres 
with energies in the range 0.1-1 GeV, there appears two 
(instead of one) enhancements of intensity, and with a 
very peculiar time profile. In the case of impact parameter 
u — 0.5, the intensity first rise relatively slowly compared 
with the immediately subsequent behavior, and suddenly 
decreases. Afterwards, this effect repeats itself in a specular 
way. This phenomenon is reminiscent of that found for 
putative matter violating some of the energy conditions 
(e.g. Safonova et al. 2001, Eiroa et al. 2001), although in 
that case, a non-zero shear is not needed in order to produce 
it and the flux rising is sharper. 

In Figure [12], we show the size and shape of the caustics 
(in the source plane) for two specific lens-source configura- 
tions and lensing parameters. We also present, for compar- 
ison, the sizes of the 7-spheres for different energies. It is 
possible to see, then, that the biggest 7-sphere here consid- 
ered, the one corresponding to E =10 GeV, is larger even 
than the caustic size for the lensing parameters used in the 
left panel of the Figure. It is for this 7-sphere that we see 
only one intensity enhancement, translated into a peak of the 
magnification. The reason is that light from different parts 
of the source smooth down the enhancement produced by 
the only part crossing the caustic. Therefore, for the biggest 
7-ray sphere, the caustic curve is seen as a point. However, 
for the smaller 7-spheres, the relative motion of the lens 
and source makes the latter to encounter two caustics, lead- 
ing to different peaks in the time profile. Contrary to the 
caustic presented in the left panel of Figure |l^, the lensing 
parameters used for the right panel produce a caustic curve 
actually formed by two separated pieces. When 8 = there 
would not be a direct crossing of the caustic, even for the 
largest 7-spheres. Two peaks profiles would be explained by 
double caustic crossing. In this case, we would see two peaks 
even for the largest 7-spheres. 



5.3 Spectral evolution 

The differential magnification makes its way to a spectral 
slope change at medium energies (see Figures ^| and [n], 
right panels). This break, that we predict as a distinctive 
feature for microlensing events of 7-ray blazars, together 
with its peculiar time evolution -it is shifted towards high 
energies as times goes by since transit- can be used to differ- 
entiate this from other phenomena. The temporal evolution 
of the break for the cases of zero shear is shown in Figure 
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Figure 10. Lensing results for a single point lens with an impact parameter u = b/r s = B/R s = 0.5 (upper panel) and 2 (lower panel). 
Left: light curves, from top to bottom E =0.1 GeV, 0.5 GeV, 1 GeV, 10 GeV. Right: spectra, from top to bottom T = vt/R E =0, 0.02, 
0.05 and 0.10. Lensing parameters were chosen as k s = 0.8 and 7 = 0. 



The spectral evolution shown by lensing with non-zero 
shear is notoriously modified. In Figure [H], for instance, the 
spectral evolution starts at T = 0, then the intensity grows 
towards higher energies because of the encounter of the right 
peak in the left hand side of the left panel plots, and then 
continues downward as the len s g ets away from the caus- 
tic. As it can be seen in Figure |ll] and successive ones, it is 
not possible to plot in this case the evolution of the spectral 
break with time as a continuous function. The break in the 
spectrum appears only once in this energy range, for the ear- 
liest time, and then presumably moves upwards, to energies 
corresponding to non-magnified 7-spheres. The evolution of 



the spectrum itself (and not of the position of the break) is 
what can help to decide in these cases if the putative obser- 
vations are associated to a microlensing event. 



6 OPTICAL DEPTH 
6.1 Time scale 

The characteristic time of a microlensing event is the time 
that the source takes to cross the Einstein radius of the lens. 
It is given by 
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Figure 11. Lcnsing results for a single point lens with an impact parameter u =0.5 (upper panel) and 2 (lower panel). Left: light curves, 
from top to bottom E =0.1 GeV, 0.5 GeV, 1 GeV, 10 GeV. Right: spectra, from top to bottom T = vt/R E =0.4, 0, 0.8 and 1.2. Lensing 
parameters were chosen as k s = 0.20 and 7 = —0.16 and 9 = 0°. 



to 




4GM DoiDi, 



Do 



(56) 



where v is the lens velocity. For a cosmological source, us- 
ing the Friedmann-Lemaitre- Robertson- Walker model of the 
universe, with z the redshift of the object, and the Hubble 
parameter Ho and the current ratio between the actual and 
the critical density of the universe flo being well constrained 
in the range 



50 



km 
s Mpc 



< H < 100 



km 

s Mpc' 



(57) 



0.1 < Q < 1, 
we have 



D(zi, Zj ) 



with 



2c (1 — Op — GiGj)(Gi — Gj 
Ho Q 2 (l + «;)(! + z 3 y 



(58) 



(59) 



Gi,j = (l + Q Zi,j) 2 - (60) 

Taking as an example, Hq = 75 km s _1 Mpc -1 , fi = 0.2, 
a lens mass of M = O.IMq, a velocity v = 5000 km s _1 , 
2 = 0.1 for the lens, and z = 1 for the source, we obtain 

to = 1.5 x 10 7 s ~ 170 days (61) 
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Figure 12. Left: Caustics for k b = 0.20 and 7 = —0.16 (diamond shaped curve). Right: Caustics for k s = 0.89 and 7 = —0.34 (triangle 
shaped curves). Sources inside the caustic have four images and those outside have two. The dashed circles are the gamma spheres, the 
bigger one corresponding to E = 10 GeV, then, 1, 0.5, and 0.1 GeV. 
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Figure 13. Temporal evolution of the spectral break in models 
with zero shear. The thinner solid lines represent linear fits, show- 
ing the general trend to increase the energy at which the spectral 
break is produced as a function of time. The impact parameter in 
this example is u = 0.5. The x-axis is T = vt/R^. 



The width of the peaks in the light curves (see Figure p"(i| ) is 
about 



At ~ 0.05t 7.5 x 10 5 s ~ 9 days. 



(62) 



These numbers would change to 1202 days and 60 days, re- 
spectively, if the lensing were produced by a 5Mq star. It is 
interesting to compare such numbers with typical timescales 
of high- latitude sources determined through the EGRET ex- 
periment. The Third EGRET catalog was constructed dur- 
ing a period of six years, dividing the total time span in 
viewing periods with a duration of about 15 days. Then, 
when the mass of the lensing object is sub-solar, each peak 
in a microlensing light curve can be completely within a 
single EGRET viewing period: the phenomenon can lead, in 
principle, to a very variable source, with 7-ray fluxes varying 
from detection to upper limits in consecutive viewing peri- 
ods. On the contrary, for a 5Mq star, a single peak would 
last several viewing periods and the 7-ray source could ap- 
pear as a steady, non-variable detection. 



6.2 Linear size of the source 

We have mentioned before that in order for a microlensing 
event to occur, the linear size of the source, x, should be less 
than the Einstein radius, 



x < 2R E 



Dos 
Dal 



(63) 



For the same typical redshifts considered in the previous 
section 



Re = 2.23 x 10 1 



M 



1/2 



(64) 
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and we can then obtain the following relationship between 
the source size and the mass of the lens: 



x < 6.37 x 10 1 



M 

M© 



1/2 



(65) 



For a source with x = 10 14 cm strong magnification occurs for 
stars with masses M/Mq ^ 2 x 10~ 3 . This makes most of the 
stars MACHO- like objects in a galaxy able to produce strong 
gravitational lensing effects upon the innermost regions of 
background active galactic nuclei. The smaller masses can 
give rise to very rapid events (Romero et al. 1995). 



6.3 Expected number of microlensing cases in 
7-ray catalogs 

The concept of optical depth, V, was originally introduced 
in gravitational microlensing studies by Ostriker and Vietri 
(1983), and it was later applied by Paczyhski (1986). If V is 
smaller than unity, it provides a measure of the probability 
of microlensing. Alternatively, F can be defined as the ratio 
of the surface mass density of microlensing matter to the 
critical mass density of the galaxy (Paczyhski 1986). The 
value of F depends on the model adopted for the matter 
distribution along the line of sight. 

It is usually assumed that the a priori probability of 
finding a small group of distant, gravitationally magnified 
objects is below 1%. Recent results (Wyithe & Turner 
2002), taking into account the clustering of stars in inter- 
posed galaxies, give for the a priori probability of finding 
magnified sources in random directions of the sky values 
between 10 -2 - 10 -3 . In those directions where there is 
gravitational lensing, the probability of having large local 
values of optical depth is high. 

The high surface mass density associated with the core 
of normal galaxies along with the usual assumption that 
most of this mass is under the form of compact objects 
naturally leads to high optical depths for microlensing. For 
instance, in the case of the lensed quasar Q2237+0305, 
where four images are well-resolved, lensing models indicate 
values of r ~ 0.5 (Schneider et al 1988; Wambsganss & 
Paczyhski 1994), which are corroborated by the detections 
of microlensing-based optical variability with relatively high 
duty cycles (e.g. Corrigan et al. 1991, Wozniak et al. 2000, 
Witt & Mao 1994). Other lensed sources display even higher 
duty cycles (e.g. Koopmans & de Bruyn 2000). 

In our case, the number of potential compact 7-ray 
emitting background sources is large. The last version of 
the Veron-Cetty & Veron's (2001) Catalog -which is still 
very incomplete at high redshifts- contains more than 

10 3 already identified blazars, in addition to more than 

10 4 quasars and other less energetic AGNs. The GLAST 
mission itself is expected to pinpoint about 10 4 7-ray 
emitting blazars with unparalleled resolution (Gehrels & 
Michelson 1999); also the number of unidentified sources at 
high latitudes is expected to be large. If the actual total 



number of 7-ray emitting blazars is in excess of, say, 10 r , (1 
blazar out of 10 000 normal galaxies) they could produce 
many of the expected detections by GFAST. Even when 
considering reduced probabilities for microlensing, scaling 
as r/A 2 with r being the local optical depth and A the 
magnification, an interesting number of detections could be 
potentially ascribed to microlensing. 

A crude estimation of the number of possible 7-ray 
sources produced by microlensing can be obtained as the 
product of three factors: random lensing probability x local 
lensing probability x number of background sources, i.e. ap- 
proximately 5 x 10 -3 x l/A 2 x 10 7 . The uncertainty in the 
previous expression, however, is large. We can only roughly 
estimate the total number of background sources, but the 
value of A they need to become visible i.e. with fluxes above 
the sensibility of EGRET and/or GLAST, will depend on 
the luminosity function of 7-ray emitting AGNs, which is 
unknown. According to the different sensitivities of both in- 
struments, we could expect, perhaps, a number of detections 
of ~ 10 and ~ 100, respectively. 



6.4 Number of events per light curves 

We shall now estimate the number of microlensing events 
expected for a 7-ray blazar with a galaxy interposed in the 
line of sight. For the galaxy mass distribution we shall adopt 
the model used by Griest (1991) -see also Eq. (^9|) above-, 
where the density profile is given by 



p{r) = Po 



a 2 + (x 2 + y 2 + z 2 ) ' 



(66) 



and where po is the mass density in the center of the galaxy 
and a is the core radius of the halo. We shall consider that 
the entire microlens population is at the same distance from 
the observer. The projected surface mass density of the mi- 
crolensing at a radius r 2 = x 2 + y 2 = constant from the 
center of the galaxy will be 

!• DC 

E(r) = / dz p 



np a 



a 2 + (x 2 + y 2 + z 2 ) 



(67) 



Assuming that E(r) is constant for the smallest impact pa- 
rameters and equal to E c = E(r = 0), we obtain 



EcTTlt R E 

M 



(68) 



with K a = E c /E cr it. Usually B is taken as unity. This would 
happen, for instance, when the distant quasar and the galaxy 
in which the lenses reside are perfectly aligned (as in the case 
of PKS 0537-441, see Romero et al. 1995 for discussion). In 
this case, all lenses are close to the center of the galaxy, and 
the previous approximation is valid. After some algebra, and 
taking derivatives, we may write this latter expression as 



dT = fiB 2 = |^^B 2 = F^ 

2_j cr it 2jc 2-*crit ^>c 



(69) 



In order to make further estimates, we need to assume some 
parameterization for the mass function of the lenses. The 
standard choice is a power law of the form N(M) tx M~ a 
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(e.g. D'Antona & Mazzitelli 1986). It is likely, however, that 
a single power law will not suffice for all different types of 
stars. The value of a, also known as the Salpeter index, is 
usually taken as 2.35 (Salpeter 1955), although more massive 
stars may require a steeper index, whereas it could probably 
be flatter for the less massive stars. As was done in Surpi et 
al. (1996), we normalize the mass function to yield a total 
surface density E c = K s S cr i t at the center of the galaxy. The 
differential surface density in objects with masses between 
M and M + dM at the center would then be 

dE c MN(M)dM 



tic 



Sm MN(M)dM 



(70) 



Using this last expression in Eq. (|f 

l- 

i / w \ 

dr = 



C(a, M n 



M 



we find 
M 



TT~ ) d — , (71) 



with 



ax ) — 



v M e ) 



in 



for a/ 2, 
U mi J ' 



for a - 



(72) 



Following Subramanian & Gopal-Krishna (1991), if the 
lenses have the same mass and the observing period is At, 
the expected number of events in the case of a background 
source moving with velocity v will be 

^_ 2 At 



N = r i + 



TlBto 

where F is the optical depth and B 
Eq. (jrij), we obtain 

2vAt 



b/R E 



(73) 
Then, using 



dN 



ttBRe 



+ 1 



C(a,M" min ,M max ) 



M 



M 



tH — (74) 



It is convenient to write the previous expression in more 
useful units. We define the dimensionless parameter 



_ 1 30days 2 Dos 
~ 7v ^/GMq- C V D ol D, 
to finally obtain 



(75) 



dN = V- 



(-) 

C 30days c \M Q J 



1/2- 



+ 



C 



M 



1/2- 



d 



(76) 



M 

For typical values, say z„ ~ 0.9 and zi ~ 0.1, T> ~ 2.1. Eq. 
(^) can now be used to obtain the number of events per 
lens mass interval as wellas the total number of microlensing 
events. 



included in the total mass range (M m i n ,M max ), during At 
days of observations is 



N 



At 

M-i—M 2 



D(Q,Afi,M 2 )F 
C 

2 16 /~~T- D ^ M ^ M ^ At 
■ VKs cC(a,M min ,M max )30days' 



(77) 



where D is given by 



D(a,M 1 ,M 2 ) = 



( M2-\ 2 a _ ( Mj\ ' 



2a 



for a # -, 



■»(t) 



for a 



(78) 



and we have assumed the previously mentioned redshifts 
(z s ~ 0.9 and zi ~ 0.1) to fix the numerical coeffi- 
cient. Clearly, the total number of events will strongly 
depend (apart from the expected influence of V and k s ) 
on the relative source-lens velocity and the Salpeter index a. 

Blazars where the bulk of the high-energy emission is 
produced in a superluminal component with apparent veloc- 
ity v > c in the lens plane will produce 7-ray sources with 
the highest levels of variability. In Figure Il4| we present three 
curves showing the total number of microlensing events in a 
5 years period for the case of the source being a blazar su- 
perluminal component, a subluminal (v — 0.5c) component, 
and for a fixed core source with a foreground lens moving 
with v ~ 1/lOOc. We find that during the EGRET observing 
time, and in the case of the lensing of a superluminal source, 
hundreds of events can be expected^] These 7-ray sources 
can strongly fluctuate from one viewing period to another, 
leading to high levels of variability. Instead, AGNs whose 
velocities in the lens plane are much smaller than c will be 
only mildly affected by lensing, yielding few events for usual 
values of T. These cases would produce 7-ray light curves 
shifting from detection to upper limits in separated viewing 
periods, even producing only one 7-ray detection in a period 
of five years. Depending ultimately on the mass of the lens, 
this unique event could produce a steady non- variable source 
(seen along consecutive viewing periods), or one which is 
seen only during a single viewing period. The peculiarities of 
each event will ultimately determine the time scales involved. 
Table |^ shows the results for the total number of events for 
different lensing parameters, whereas Table |^ shows the dis- 
tribution of the events by lensing mass. It is interesting to 
note that most of the events are produced by MACHO-like, 
sub-stellar mass objects. These objects present small Eintein 
radius, indicating that only the innermost 7-spheres will be 
magnified. In addition, this would entail smaller time scales, 
and within the galactic model here considered, more variable 
7-ray sources. 



For the assumed power law mass distribution of lenses 
between M m i n < M < M max , the number of expected mi- 
crolensing events by stars with masses in the range (Mi, M2), 



1 "Events" refers to single spikes within the light curves whereas 
"sources" to the number of EGRET detections in different parts 
of the sky that could be ascribed to this mechanism. 
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Table 2. Total number of events for different lensing parameters 
and Salpeter indices produced in a period of five years. Two cases 
(bottom panels) show the results for very low optical depths: even 
in that cases, significant number of events are expected if the 
sources have apparently superluminal velocities. 
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Figure 14. Number of events as a function of the Salpeter index 
a for the case of apparently super and sub-luminal movements 
of the source in the lens plane. Lensing parameters are 7 = 0, 
k s = 0.4, and T = 0.2. For v -C c, the time scale of one event 
is comparable with EGRET observing time. For v ~ c the time 
scale of a single event is comparable with one EGRET viewing 
period. 



Table 3. Distribution of the number of lensing events by mass 
for a galaxy with n s = 0.4, 7 = 0, and T = 2. The Salpeter index 
is a = 2.1 for the first panel, a = 2.3 for the second, and a = 2.5 
for the third one. The relative source velocity in the lens plane is 
chosen as v = 3.5c, 0.1c and 0.01c, from top to bottom. 
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Figure 15. Light rays are traced back towards the source. Each 
of them is deflected in the lens plane. Different parts of the lens 
plane can generate deflections converging onto the same pixel in 
the source plane. There, the magnification map -the number of 
rays per pixel- is made. Once this is done, a source with a given 
size moves in this magnification plane and generates a light curve. 



6.5 Critical assessment 

Although the Chang-Refsdal scheme used in the previous 
section is a very important tool in the understanding of the 
model herein proposed, it cannot be used to claim precise 
predictions for the magnifications. For high values of k s and 
7 a more complex caustic pattern should be considered in 
order to get reliable light curves and the correct magnifica- 
tions. The quantitative predictions for the amount of mag- 
nification and the form of the light curves could actually 
change in more detailed simulations of a galaxy core (where 
the effects of thousands of stars are considered simultane- 
ously). The rest of this paper is devoted to develop full nu- 
merical models of the situation presented from an analytical 
point of view in the previous sections . 



7 MAGNIFICATION MAPS 

As we have seen, the parameters that describe a microlens- 
ing scenario are the dimensionless surface mass density k 
-expressed in units of the critical surface mass density- 
and the external shear 7 (cf. Kayser, Refsdal & Stabell 



1986; Schneider & Weiss 1987). The former -often also 
called convergence or optical depth- describes the amount 
of matter in front of the source. The latter is a tidal force 
caused by matter outside the light bundle. In order to 
simulate the effect of a particular combination of k and 
7, point lenses are distributed randomly according to the 
given surface mass density. If we would replace each point 
lens by a disk with radius equal to its Einstein radius, the 
total fraction of the sky that is covered by the sum of all 
these disks is equal to the optical depth. 

For the ray-shooting simulations, a large number of light 
rays (of order 10 9 ) are followed backwards from the observer 
through the field of point lenses. The number of rays de- 
termine the resolution of the numerical simulation. These 
rays start out in a regular (angular) grid. In the lens plane, 
the deflections according to the individual lenses are then 
superposed for each individual ray i, 



3=1 



AG 

„2 



3=1 



(79) 



Here Mj is the mass of the point lens j, is the projected 
vector distance between the position of the light ray i 
and the lens j, and m is its absolute value. It is the 
computations of these individual deflections what requires 
most of the time of the simulation. The effect of the 
external shear is included as well. The deflected rays are 
then followed further to the source plane. There, they are 
collected in small pixels. The number of rays per pixel (on 
average ~ 100 for a region typically of 2500 x 2500 pixels) 
is proportional to the magnification at this position. A 
two-dimensional map of the ray density -a magnification 
pattern, also referred to as caustic pattern- can then be 
produced. The magnification as a function of (source) 
position is indicated by colors. Sharp lines correspond 
to locations of very high magnification, i.e. the caustics. 
In Figure |l5| we schematically show the geometry of the 
shooting technique. 

The minimal number of lenses that have to be consid- 
ered depends on the focusing and shear values, as well as 
on the ratio between the diffuse and the total flux. The dif- 
fuse flux (e) is that coming from rays that are deflected into 
the receiving area from stars far outside the region where 
microlenses are considered, and should be consistently low. 
An approximated expression for the number of lenses to be 
included in each magnification map is (Wambsganss 1999) 



N„ 



3f^ 



(1-k) 2 -7 2 e 



(80) 



which entails values from several hundreds (for n < 0.4) 
up to several hundred thousands (for k ~ 1) stars, in the 
case of zero shear and e = 0.01. Since for computing the 
map with high resolution, this would imply literally several 
million billions of operations, an intelligent numerical 
routine should come to the rescue to make this problem 
feasible. The technique is called hierarchical tree code. 
Basically, it groups the stars as a function of their distance 
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to the light ray. Since their influence decreases as gravity 
does, with r~ 2 , one does not need to compute the deflection 
from all the stars involved with the same level of precision 
to obtain an overall sensitive result. A detailed account 
of the numerical technique can be found in the paper by 
Wambsganss (1999). 

In general, a smooth out distributed surface mass den- 
sity k c also contributes to the deflection as well, and the 
general microlensing equation to be solved is 

y= ( 1_ o~ 7 i-° + 7 (si) 

In observational situations n c and 7 are obtained from 
macrolensing simulations of the resolved source. Rays repre- 
senting a square (x\ = X2) are then mapped onto a rectangle 
with a side ratio T = (1 — k — K c — y)/(1 — K — k c + 7). But, 
as we would like the receiving -and not the shooting- area 
to be a square (the pixel) the shooting field (i.e. the area in 
the lens plane in which rays are mapped) is chosen to be a 
rectangle of size T _1 . 

The fact that source, lens and, observer are moving 
relative to each other, means that the source is affected 
by a variable magnification as a function of time: it moves 
through the magnification pattern and so we measure a 
variable flux. When a source crosses a caustic, formally 
two very bright new (micro-)images appear or disappear. 
However, their angular separation is much smaller than the 
resolution of any telescope. Therefore, we can only measure 
the combined total brightness, which can produce dramatic 
'jumps' in the observed flux. Assuming geometrical optics 
and a point source, the amount of the magnification 
is infinite. However, since any real source is finite, the 
magnification stays finite as well. The exact amount of 
magnification depends on the source size: the smaller it 
is, the higher the magnification. In order to determine the 
light curve for a finite source, its brightness profile has to be 
convolved with this two-dimensional magnification pattern. 



7.1 Results of the numerical simulations 

In Figure [l6] we show the magnification map for the case 
in which k = 0.5 and 7 = 0.0; the brighter the region, the 
stronger the magnification. The characteristic critical lines 
of the Chang-Refsdal model appear in this map, as well as 
the diamond shaped structures due to the close-by star. 
In the bottom left of the panel, we show the size of the 
source for three different energies. The innermost pixel (the 
center of the circles) represent the size of the lowest energy 
7-sphere, that corresponds to E = 100 MeV. The second 
circle is the size of the 1 GeV 7-sphere, whereas the largest 
one is the corresponding size of the E = 10 GeV emitting 
region. It is clear the origin of the differential magnification: 
the 7-spheres will be affected according to their size while 
moving in the caustic pattern. 

The numbered lines in Figure hq represent different 
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Figure 18. Light curve for trajectory number 7, producing one 
of the highest levels of magnification out of map given in Figure 
The line coding is as in Figure [l?]. 

source trajectories. Numbers 1 to 6 are trajectories common 
to all maps we shall present, and are defined a priori in 
our code. The light curves for each of these six trajectories 
are given in Figure ^ We can see that the overall effect 
described with the Chang-Refsdal model appears here as 
well. In lighter color we show, for each trajectory, the 
corresponding light curve for the E = 100 MeV emitting 
region. Darker lines correspond to regions emitting photons 
of 1 and 10 GeV. The latter, in all cases, are smoother 
versions of the former and, always, the magnification is 
weaker for these regions. There is a typical factor of 10 
more magnification for the innermost regions than for 
the larger 7-spheres. We see that for these values of k 
and 7 it is possible to get a typical enhancement of 10 
times the unlensed intensity at 100 MeV. The drop to 
zero at the corners of the diagram for all light curves 
corresponding to the largest 7-sphere is an artifact of 
the simulations: the code assigns zero magnification when 
more than half of the source is out of the magnification map. 

The shape of the light curves is also worthy to comment. 
As an example, we take trajectory 3. It starts in a region of 
low magnification and continues upwards, crossing a region 
of relatively high enhancement, where for the innermost 
7-sphere the magnification is 9 times the intensity of the 
unlensed source. There are four caustic crossings there. 
We can see that for the innermost 7-sphere, the caustic 
crossings are well separated events, so we have four peaks 
in the light curve corresponding to E — 100 MeV. However, 
for the larger 7-spheres the peaks are smoothed down. We 
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Figure 16. Magnification map for lcnsing with parameters re = 0.5 and 7 = 0.0. For details, see text. The corresponding file is 
magpat k 50 g 00.gif. 
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Figure 17. Light curves for different source trajectories. Numbers corresponds to those given in Figure [^| Darker lines corresponds, 
respectively, to regions emitting photons of 100 MeV, 1, and 10 GeV, and whose emitting sizes are depicted in the bottom left corner of 
Figure ha, being the innermost point the less energetic 7-ray sphere. 



see only two broadened peaks for E — 1 GeV, and only one, 
with almost nil magnification, in the case of E = 10 GeV. If 
the 7-ray AGN is within the observing sensitivity, we would 
see a distinctive effect during the microlensing event due to 
the different sizes of the source at the different energies. 

The x-axis in Figure [I?] is a linear length scale, the Ein- 
stein radius of a solar mass star, Re(Mq) — 2.23 x 10 16 
cm. It can be translated into a time scale as t — Re(Mq)/v, 



where v is the relative velocity of the source with respect to 
the lens, projected onto the source plane. Typical time scales 
for microlensing events where discussed in Section 6.1. To 
give an example of the similar time scales predicted with a 
full caustic pattern plot, we write the time scale as 

t= R E (M e)= 0_m3 

v v/c J v ' 

With a relative velocity v — 5000 km s _1 , the length of 
each axis in Figure llq is equivalent to 14 years. The other 
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important time scale involved in a microlensing event is the 
rise time to a peak of maximum magnification. This will 
depend on the size of the different 7-spheres, and is given by 

R 0.023 R 

r = 7 = ^7MMo) yr ' (83) 

where, again, we have scaled it with the Einstein radius 
corresponding with one solar mass. Then, the innermost 
7-spheres (having R/ Re ~ 1/100) will have a rise time scale 
of about 5 days, well within an observing EGRET viewing 
period. The largest 7-spheres, with R ~ Re, can have a rise 
time of about 1 yr. Higher (lower) velocities would imply 
lower (higher) time scales. 

In order to explore the maximum possible magnification 
that this caustic pattern produces, we have selected an 
extra trajectory that crosses exactly over a conjunction 
of several caustics (trajectory 7 in Figure |l6| ). The light 
curve is separately shown in Figure hq. We can see that 
magnifications up to 45 times the unlensed intensity of 
the source are possible in this configuration. The effect of 
the differential enhancement is notoriously clear for this case. 

In Figure [19] we show the magnification map corre- 
sponding to a higher value of focusing, k = 0.8. Still in this 
case, the shear 7 is taken equal to zero. As it was found 
by Schneider & Weiss (1987), the critical structure become 
more complex with increasing k and is no longer possible 
to identify a constellation of compact objects. In addition, 
we also see (as in Figure |l6|) the tendency of the caustic 
structure to cluster, generating some crowed critical regions 
and some others devoid of high magnification patterns. The 
explanation for this was already given by Schneider & Weiss 
(1987): The clustering of caustics is just the non-linear 
enhancement of random (Poisson) clustering of the positions 
of lenses in the source plane. Over-dense regions tend to 
attract other over-dense regions, because gravity is a long 
range force. 

The most important feature shown in the map of 
Figure |l^ is that the magnitude of the magnification has 
been typically reduced with respect to the k = 0.5 case, and 
this reduction reach two orders of magnitude with respect 
to the single events. The density of caustics is so large that 
the light curve is continuously affected by them, producing 
a less dramatic combined effect. This effect was first studied 
by Deguchi & Watson (1987): The total magnification is 
always high, but the fluctuations decreases beyond k = 0.5. 
This effect can be seen in the light curves presented in 
Figure |5(]. Even for the particularly chosen trajectory, 
number 7, we find a maximum magnification of 7. The full 
light curve is presented in Figure Ell Again, the differential 
effect is notorious. 

In Figure [22] we show the magnification pattern for the 
case k = 0.2 and 7 = 0.16. The presence of shear modifies 
qualitative the magnification map. Here, most of the map 
is devoid of magnification, and so, many of the common 
trajectories (numbers 1 to 6) cross large regions of very low 
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Figure 21. Light curve for trajectory number 7, producing one 
of the highest levels of magnification out of map given in Figure 
|xs| . The line coding is as in Figure [l?]. 

or nil magnification (see particularly trajectory 6). However, 
those trajectories actually crossing the caustics produce 
enhancements in intensity typically between 10 an 20 times 
the unlensed value. These effects can be seen in the six 
panels of Figure where the differential magnification for 
the different 7-spheres is also noticeable. In addition, the 
enhancements of intensity are usually well separated (see 
for instance trajectory 4). The enhancements themselves are 
reminiscent of those produced by single events considered 
in Section 6, in those cases where the presence of shear was 
significant. 

From the point of view of unidentified 7-ray sources, 
interposed galaxies with low values of k and 7 are probably 
the most interesting case for the application of the model. 
In Figure [24] we show the case for the hand-selected light 
curve. In that case, the magnification reaches a factor of 
65. The probability for this trajectory is less than those of 
average enhancement. However, even if the probability is 
reduced by, say, a factor I /A 2 , it is possible to expect many 
cases of high magnification, like that presented in Figure 

H 

Finally, in Figure we show the magnification map 
corresponding to the case K = 0.9 and 7 = 0.4. Again, the 
high value of k makes the critical structure highly complex. 
Typically, the magnification values are below a factor of 10 of 
the unlensed intensity, although some trajectories are found 
(see, for example, curves number 3 and 4 of Figure |2(j) where 
a factor of ~ 10 is reached in two well separated regions. 
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Figure 19. Magnification map for lcnsing with parameters k = 0.8 and 7 = 0.0. For details, see text. The corresponding file is 
magpat fc 80 9 00 . gi f . 
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Figure 20. Light curves for different source trajectories. Numbers corresponds to those given in Figure |l9| The line coding is as in Figure 
0- 



Trajectory 7 (whose light curve is given in Figure |l^) shows 
an enhancement of 22 times the unlensed intensity. 



7.2 Speculations on astrophysical applications 

Let us assume that a correct alignment of source, lens, and 
observer has been produced and that microlensing is oper- 
ative. A possible astrophysical application of the effects we 
have discussed in the previous sections is to constrain the 
exponent in the relationship between size and energy for 7- 



ray spheres in the AGNs , R oc E a . The smaller the source, 
the higher is the peak 7-ray luminosity, and shorter will re- 
sult the rise time scale (r = R/v). Therefore, by observing 
throughout the peak at different energies, one could deter- 
mine the relative size of the source at such energies, and then 
test the radius-energy theoretical relationship. This could be 
done within the range of GLAST capabilities. If, somehow, 
we know the relative velocity and the redshifts, we could in 
even determine the sizes of the different emitting regions. 
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Figure 22. Magnification map for lensing with parameters k = 0.2 and 7 = 0.16. For details, see text. The corresponding file is 
magpat k 20 g 16.gif. 
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Figure 23. Light curves for different source trajectories. Numbers corresponds to those given in Figure |22| The line coding is as in Figure 
0- 



This would directly impact on the underlying 7-ray models 
(e.g. Becker & Kafatos 1995, Blandford & Levinson 1995). 



8 SUMMARY AND FINAL COMMENTS 

To summarize, we have shown in this paper that 

• Some of the high-latitude unidentified 7-ray sources 
(both variable and non-variable) could be weak 7-ray 
emitting AGNs that are magnified through gravitational 



microlensing by stars in foreground galaxies. 

• Although small, the probability of gravitational mi- 
crolensing could be enough to explain a handful of the 
EGRET detections, and maybe many of the forthcoming 
GLAST detections. 

• During a 7-ray variability event produced by microlens- 
ing of a blazar, there is a peculiar spectral evolution that 
could be detected in principle by the next generation of 7- 
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Figure 25. Magnification map for lcnsing with parameters re 
magpatj.90g40.cji/. 



0.9 and 7 = 0.4. For details, see text. The corresponding file is 
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Figure 26. Light curves for different source trajectories. Numbers corresponds to those given in Figure |25| The line coding is as in Figure 
0- 



ray observatories with fine spectral capabilities, like GLAST, 
or even by ground 7-ray telescopes located at sufficient alti- 
tude, like 5@5 (e.g. Aharonian et al. 2001). 

It should be clearly stated that this model can not account 
for all unidentified 7-ray detections at high latitudes. 
However, it is interesting to ask whether the proposed mi- 
crolensing scenario could be responsible for 7-ray variability 
of some radio loud AGNs already detected by EGRET. In- 
deed, there is one possible case, related with the source 3EG 
J1832-2110, which has been identified with PKS 1830-211 



(Mattox et al. 1997, Combi & Romero 1998). The latter is 
a flat-spectrum radio source, proposed to be a gravitational 
lensed QSO by Pramesh Rao & Subrahmanyan (1988). The 
7-ray source is likely variable, presenting a value of / = 2.5, 
and a steep spectral index, T = 2.59 ± 0.13. Both facts, 
variability and a steep spectra, argue against a galactic 
origin for this source. Mattox et al. (2001) assign to this 
pair an a priori probability of 0.998 of being correct, and list 
it within the most likely AGN identification of the Third 



© 0000 RAS, MNRAS 000, 000-000 



28 Torres, Romero, Eiroa, Wambsganss, & Pessah 



65 
60 
55 
50 
45 

S 40 

o 

S 35 

fi 
tic 

E 30 

O " 

c 

20 
15 
10 
5 





2 4 6 8 10 

length scale (unit: Einstein radius of solar mass star) 



Figure 24. Light curve for trajectory number 7, producing one 
of the highest levels of magnification out of map given in Figure 
p^ . The line coding is as in Figure [l?]. 



EGRET Catalog. 

High resolution radio images obtained from several 
interferometric arrays have revealed that the source has 
a ring-like structure with two bright components on sub- 
arcsecond scales (Jauncey et al. 1991). This suggests a close 
alignment of the lensed source behind the lensing object. 
Two absorption systems have been detected at z ~ 0.89 
(Wiklind & Combes 1996) and z ~ 0.193 (Lovell et al. 
1996), so it seems likely that the image of the background 
QSO (with a redshift z > 0.885, Mattox et al. 2001) is 
lensed by two different extragalactic objects, what would 
undoubtedly enhance the optical depth and the number of 
expected single microlensing events. 

Indeed, Combi & Romero (1998) have already proposed 
that the 7-ray emission of 3EG J1832-2110 (then 2EG 
J1834-2138) could be produced by gravitational microlens- 
ing, using exactly the same ideas we have deepen in this 
paper. They have found that assuming a redshift z 3 — 1 
for the background source and Zi ~ 0.89 for the lens, a 
MACHO-like object with M ~ 0.02M Q and moving with 
a low velocity of only v ~ 1000 km s" 1 would be enough 
to produce the observed variability. For this to be possible, 
the size of the 7-ray emitting region should be of about 1.5 
10 15 cm, in good agreement with the source sizes used in 
this paper. These results can be slightly modified by new 
measurements of the quasar redshift, but will not change 
substantially (Oshima et al. 2001). 
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Figure 27. Light curve for trajectory number 7, producing one 
of the highest levels of magnification out of map given in Figure 
The line coding is as in Figure [nj. 



It is likely then that the first realization of this pro- 
posed mechanism have been already observed for the pair 
3EG J1832-2110/PKS 1830-211. One thing should be re- 
marked, though: in this case the background source is al- 
ready a strong radio emitter -what indeed facilitates the 
identification-. This was not the general case we have consid- 
ered here, where the sources are weak enough as to yield no 
significant lower energy counterparts. A complete microlens- 
ing model for the 7-ray variability of PKS 1830-211, based 
on the discussion presented in this paper, will be presented 
in a forthcoming publication. 
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APPENDIX: SOLUTIONS TO THE ID LENS 
EQUATION 

Replacing Eqs. @ and @ in Eq. @ we obtain 
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Let us now analyze different situations, starting with the 
case e = 1, Y > 0. For it we get, 
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1 / (F + vTHi) 4 
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If we have the case e = 1, Y < 0, we get 
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For the case in which e = —1, Y > 2, the magnification 
results in 
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And lastly, for the case in which e = —1, Y < —2, 
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The four cases can be jointly written as 
1 
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16- (eY - VY 2 +4e) 4 / 
Simplifying the expression between the parentheses, we have 
(eY + v / Y 2 + 4e)' ; 
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Then (using that e 2 = 1), we obtain 



1 Y 2 + 2e 

sff(Y) 
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|1 - «| s 5 (Y)YVY 2 + 4e' 
and as sg(Y)Y = |Y|, Eq. (§|) follows. 
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